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Abstract 

The cross section for the elastic photoproduction of mesons (7^ p^p) has 
been measured with the HI detector at HERA for two average photon-proton centre- 
of-mass energies of 55 and 187 GeV. The lower energy point was measured by ob- 
serving directly the p^ decay giving a cross section of 9.1ib0.9 (stat.)ib2.5 (syst.) ph. 
The logarithmic slope parameter of the differential cross section, dfi/dt, is found to 
be 10.9 lb 2.4 (stat.) it 1.1 (syst.) GeV~^. The p° decay polar angular distribution is 
found to be consistent with s-channel helicity conservation. The higher energy cross 
section was determined from analysis of the lower part of the hadronic invariant 
mass spectrum of diffractive photoproduction and found to be 13.6 it 0.8 (stat.) ± 
2.4 (syst.) ph. 
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1 Introduction 

The photon, exhibiting both point-hke and hadronic behaviour, has a rich structure that is 
still not yet fully explored. New opportunities for studying photon interactions are opened 
by the HERA electron-proton collider. The nearly real photons produced by HERA's 
electron beam extend the energy range for photoproduction processes by more than one 
order of magnitude, up to the equivalent of tens of TeV in fixed target experiments, 
providing a large kinematic range for studying the properties of the photon. The HI and 
ZEUS collaborations have begun to explore this domain, e.g. with measurements of the 
total 7p cross section [||, 0, ^ and vector meson production 0, ^, ||, ||, 0]- 
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In this paper we report the HI results for the elastic p° photoproduction cross sec- 
tion at two average total energies, (W^p) = 55 and 187 GeV, of the jp centre-of-mass 
system (cms). For the low energy measurement the associated differential cross sections, 
dcr/dM^^, where M^^^^ is the invariant mass of the p° decay products, da/dt, where t is 
the square of the four-momentum transfer at the proton vertex, and the polar decay angle 
distribution, dcT/dcos6'*, of the p° are also presented. 

At lower energies the photon is known to exhibit behaviour resembling that of a 
hadron. In particular, the reaction 

7 + P^P°+p (1) 

at cms energies of 10-20 GeV behaves much like elastic scattering of hadrons: the cross 
section is roughly constant with cms energy, the dependence on t of da/dt is close to 
exponential and s-channel helicity conservation holds (the retains the helicity of the 
photon). This is interpreted in the Vector Meson Dominance (VMD) model [§, Q as the 
fluctuation of the photon into a virtual p° meson, with subsequent elastic scattering of the 
p° on the proton. Regge theory |TD| inspired models have been used successfully to describe 
the energy variation of the total cross section, with a 'soft' pomeron term dominating at 
high energies |TT|. Models based on perturbative quantum chromodynamics (pQCD) may 



combine a soft VMD component with perturbatively calculable components to arrive at a 
total 7p cross section. The 7p — > p^p cross section may then be completely determined by 
the parameters used to fit the total jp cross section. Thus studies of reaction (^ at HERA 
energies help to constrain such models, as well as to clarify the nature of the pomeron. 
Thorough exploration of these interactions in the HERA energy regime is important for 
an eventual understanding of the hadronic properties of the photon. 

Almost real photons, radiated at small angles by HERA's electron beam (26.7 GeV 
electrons in 1993 and 27.5 GeV positrons[| in 1994), can interact with the 820 GeV protons 
at very high energies in the jp cms. In the VMD model, one expects dominant production 
of the vector meson series u, along with other hadronic production. Elastic p^p 
interactions are characterised by a scattering of the virtual p° meson and the proton at 
low \t\, leading to a p° at low pt relative to the HERA beam axis. The p° subsequently 
decays into a tt+tt^ pair. 

The kinematic variables are defined as follows: 

t = [q-vf = [p-p'f 

q ■ P 
y = kTp 
s = {P + kf 

where y is the inelasticity parameter, is the negative square of the four-momentum 
transfer at the electron vertex, k{k') denotes the four vector of the incoming (scattered) 
electron, P{P') of the incoming (scattered) proton, q of the photon and v of the vector 
meson. 



^The term 'electron' will subsequently be used to refer to both electrons and positrons. 



4 



The two measurements presented in this paper were made under very different ex- 
perimental conditions. In the low W-yp analysis, where the p° rest frame is boosted by 
about 1 unit of rapidity toward the electron beam direction, the p° decay products were 
observed in the central tracking detector (see section ^ for details of the HI detector), 
enabling good reconstruction of the hadronic final state (excluding the proton, unseen in 
both analyses). The W^p range (40 < W^p < 80 GeV) is, however, too low for the scat- 
tered electron to be detected in the small angle electron tagger. The kinematic quantities 
were therefore reconstructed from the p°'s two decay pions, using the assumptions of an 
elastic interaction and of = GeV^. The resonant contribution to the cross section is 
deduced from the analysis of the line shape. 

The high W-yp analysis used those data with an electron tag so that the W,yp range 
was 164 < W-yp < 212 GeV. At these high values of W-^p the p° rest frame is boosted 
by about 3 units of rapidity toward the electron beam direction and the decay products 
are kinematically forced outside the acceptance region of the tracking detectors, requiring 
reconstruction of the hadronic final state in the calorimeters. Since the spatial and energy 
resolutions of the calorimeters are worse than those of the trackers no measurement was 
possible of the resonance line shape, t dependence, or polar decay angle distribution. 

The measured ep cross section is converted into a jp cross section, averaged over the 
available W.yp and Q"^ range so that (cr^p) = / ^^^q^ dydQ^/ J $(?/, Q^) dydQ^. In order 
to evaluate the photon flux, ^{y,Q^), we ignore any lepton beam polarisation and the 
longitudinal (helicity zero) contribution is neglected. In this case the flux for transverse 
virtual photons is given by the Weizsacker- Williams approximation which we take in the 
form: 

^^r>2\_ ^ ( l + jl-yr 2(1 -y) Ql\ . ^2 .2) 



where a is the fine structure constant and mg is the electron mass |jT2|, |1^. The quantity 
Qmin is the minimum photon virtuality which is kinematically allowed. The upper limit on 
the range over which the measurements were taken, Q^ax' determined by the experi- 
mental procedure used. For the low W-yp analysis the requirement of |trec| < 0.5 GeV^ (see 
section |) implies effectively Ql^^^ ^ 0.5 GeV^ and the mean is (Q^) = 0.035 GeV^. 
The requirement of the scattered electron to be detected in the tagger in the high W.yp 
analysis gives Q^^^ ^ 0.01 GeV^ and (Q^) = o.OOl GeVl 

With an assumption for the dependence of the virtual photon-proton cross section, 
the total elastic p° cross section at = may be extrapolated from (cr-yp). In a VMD 
model in which the p° dominates, the longitudinal and transverse virtual photon- 
proton cross sections are related by 0"^*^ = a^*p ■ Q^/M"^ and evolve such that a'^*p oc 
1/(1 + Q'^/MpY. With these assumptions, corrections averaged over the selected samples 
of events amount to +3.9 % and +0.2 % for the low and high W-yp analyses respectively. 
All results quoted in this paper do not include any correction for this extrapolation. 

2 The HI Detector 

A complete description of the detector is available elsewhere . Only a brief description 
of the parts of the detector relevant to these analyses is included here. In HI the z-axis 
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points in the proton beam direction, denoted by the term 'forward'. 

The HI Central Tracking Detector (CTD) provides charged track measurement and 
triggering in the pseudorapidity [r] = — lntan^/2) range —1.5 < 7] < 1.5. Its main 
components are two concentric drift chambers, Central Jet Chambers 1 and 2 (CJCl 
and CJC2). Located at the inner radius of CJCl and CJC2 are two sets of cylindrical 
drift chambers for measurement of the 2;-coordinate and two proportional chambers, CIP 
and COP. The entire CTD is situated within a 1.15 T solenoidal magnetic field, parallel 
to the proton beam axis. The transverse momentum resolution for charged tracks is 
cr{pt)/pt ~ 0.009 ■ pt{GeV) © 0.015. The CTD was the detector used to trigger and 
reconstruct the two particle hadronic final state in the low W^p analysis. 

The CTD is surrounded by a finely segmented Liquid Argon calorimeter (LAr) that 
provides hadronic and electromagnetic energy measurements in the range —1.5 < rj < 
3.3. Instrumented sections of the return yoke of the magnet outside the LAr calorimeter 
provide tracking and hadronic 'tail-catching' calorimetry in the range —2.5 < rj < 3.4. 
These detectors were used to veto inelastic and cosmic muon events. 

In the forward region two detectors were used to veto proton dissociation events: 
the forward muon detector (FMD) and the proton remnant tagger (PRT). They were 
used in the high W^p analysis to tag hadrons produced at large values of pseudorapidity 
(5.0 < r] < 7.5) by detecting secondary particles issued from collisions with the beam pipe 
or adjacent material. The FMD consists of planes of drift chambers at 6 < 2; < 9 m and 
the PRT is made up of scintillator planes at z = 24 m. 

Two detectors provide coverage in the backward direction. The Backward Electro- 
magnetic Calorimeter (BEMC), a lead scintillator calorimeter covering the range —3.4 < 
7] < —1.5, was used to reconstruct the hadronic final state in the high W^p analysis, while 
the Backward Proportional Chamber (BPC), —3.0 < 1] < —1.5, helped in the rejection of 
beam-gas background. 

Behind the BEMC are located the Time of Flight (ToF) scintillator walls, covering a 
pseudorapidity range of —3.5 < r] < —2.0. Their timing information was used to reject 
background upstream proton interactions and in the high W-^p analysis as a trigger. 

The luminosity system, monitoring the reaction ep —>■ e'-fp, consists of two TlCl/TlBr 
crystal calorimeters. The electron tagger is located ai z = —33 m and the photon detector 
at z = —103 m. The electron tagger was also used to detect and trigger on the scattered 
electron for the high W^p analysis. 



3 Monte Carlo 



Two Monte Carlo (MC) generators were used in these analyses, PYTHIA |T5| and PHO- 
JET |jl6|. Five MC photoproduction mechanisms are distinguished: elastic interactions 
(EL) which describe 7p Vp where V stands for one of the vector mesons p°, u or 0, sin- 
gle photon diffractive dissociation (CD) in which the photon dissociates (7^ Xp), single 
proton diffractive dissociation (PD) in which the proton dissociates (7^ VY), double 
diffractive dissociation (DD) in which both photon and proton dissociate {'jp XY) and 
non-diffractive interactions (ND). This last category includes all reactions not included 
in the previous classes. The ratios of the cross sections for the 5 processes were chosen 
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following comparisons made using the high W^p data (see section |^) and the constraint 
imposed by the measurement made of the total W^p cross section 0: 

Cel '■ ctgd '■ CTpB : (Tdd : ctnd = 1 '■ 1-25 : 0.75 : 2.00 : 5.0 . (3) 

All channels were generated with a dependence following equation ^ with no de- 
pendence on W^p for the photon-proton cross section. 

Both MC generators have similar inputs for the EL channel. Events are generated with 
a t-dependence of e** where 6 11 GeV~^, a decay angular distribution in the p° rest 
frame cos^* proportional to sin^^*, as expected for a transversely polarised (helicity 
±1), and a skewed mass distribution consistent with that measured (see section P). The 
relative rates of p^, uj and cf) production were chosen to be 13:1.5:1 following measurements 
made at lower energies p. 

For the singly diffractive dissociation channels (GD, PD) PYTHIA assumes a t- 
dependence of 6 ~ 5 GeV^^. PHOJET assumes a h parameter that depends on the 
mass of the diffractive system, Mx-, so that 6 ^ 11 GeV~^ at Mx = Mp and becomes 
smaller for higher Mx |1^. For the double diffractive dissociation and for both MC 
generators, h — ?> GeV~^. 

Diffractive dissociation of the proton or photon is described in PYTHIA according to 
a distribution such that roughly dN/dMx ~ l/^^x with an enhancement at lower masses. 
Masses are generated starting at 0.2 GeV (twice the pion mass) above the mass, Mj„, of 
the incoming particle, with the p° mass used for the incoming 7. PHOJET in contrast 
generates a mass spectrum according to dA^/dM|^ ~ 1/(M|^ — Mf^) at higher masses and 
treats the lower part as a single resonance. The spectrum starts at twice the pion mass 
above Mj„. 

The MC event samples included a full simulation of the HI detector and were subject 
to the same reconstruction as each of the two data samples. 



4 The Low W^p Analysis 

The data for this analysis were taken in 1993 when HERA operated with 90 bunches of 
820 GeV protons and 94 bunches of 26.7 GeV electrons with 84 colliding bunch pairs. A 
small number of non-colliding electron and proton bunches, called 'pilot' bunches, were 
included in order to estimate the beam-gas background. The data sample was limited 
to the largest continuous segment of runs containing stable and uniform central tracker 
conditions and amounted to a total integrated luminosity of 19.8 ± 1.0 nb~^. 

The most favourable trigger capable of distinguishing the soft, low multiplicity hadronic 
final state from the high level of background was the Drift Chamber r-0 trigger (DCR0) 
which distinguished individual tracks by matching CJCl and CJC2 wire hits with pre- 



defined masks The major drawback to this trigger was that its efficiency threshold, 
at pt ~ 450 MeV, was located above the peak of the decay pion pt distribution. It was 
found that the best trigger condition required one and only one track pattern identifica- 
tion by the DCR0. A requirement of two (or more) tracks above the threshold would have 
severely reduced the trigger efficiency and resulted in an unacceptable rate of background 
triggers. Thus, the second track in the event would usually have a pt below threshold and 
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would not be detected by the trigger. In order to reduce the rate of beam-gas triggers, the 
overall trigger required in addition a single vertex reconstruction by CIP and COP |jl9[ {z- 
vertex) and was vetoed if there was significant energy present in the forward part of LAr, 
excessive activity in the backward portion of CIP or an out-of-time ToF signal indicating 
a background event. 

Approximately 1.4 x 10^ events satisfied the trigger. Analysis of the data from the 
pilot bunches indicated that the vast majority of triggered events consisted of proton 
beam induced background. The data were reduced in a three stage procedure consisting 
of computer selections and a visual scan. 

Since this analysis relies on the identification of the two particle final state, it is 
particularly sensitive to the quality of track reconstruction. Having observed differences 
between the data and MC track multiplicities, a visual scan early in the analysis was 
found to be necessary in order to establish equivalent data and MC samples of clean two 
particle events. This procedure also allowed a measurement of the event reconstruction 
efficiency independent of the MC simulation. 

The first stage in data reduction employed a loose track requirement which was es- 
sentially independent of track reconstruction quality. Its purpose was to eliminate the 
obvious events that were not composed of only two particles, thereby reducing the data 
sample to a manageable size for scanning. Those events with a total of 1 track and no 
other track segments, 2 like sign, or more than 3 well reconstructed tracks were rejected. 
All other events were kept. A sub-sample of 600 of the rejected events were scanned and 
only one event was incorrectly rejected due to an improper track reconstruction (a 0.3 % 
error). 

The reduced sample of ~ 47 x 10^ events was then scanned using the HI event display 
program. Each event was required to have two opposite sign track patterns pointing to a 
common vertex and extending radially out to at least CJC2 (r > 43 cm). Only digitised 
hit information from the CJC detectors was used in the scan. Obvious cosmic events with 
tracks in the muon detector were removed. A total of 12,200 valid two-particle, opposite- 
sign events passed the scan including 623 from the proton pilot bunch crossings and 13 
from the electron pilot bunch crossings. 

Using the track fitting of the standard HI reconstruction program, the third stage of 
the reduction focussed on filtering two-particle events with well-reconstructed vertices. 
The efficiencies for data and MC events to pass this stage in the analysis were found to 
be 0.8430 ± 0.0034 (stat.) and 0.8807 ± 0.0022 (stat.), respectively. The ratio of the data 
to MC efficiencies was used as a correction factor to account for this remaining difference 
in event reconstruction (C in equation |). 

An accurate reconstruction of kinematic variables was possible using only the recon- 
structed p° and the assumption of an elastic event with Q"^ = GeV^. Figures |l|a and |l|b 
compare the reconstructed to the generated values for W-yp and t. The mean differences 
were not significant: AW^p ^ 0.4 GeV with an RMS of 1.2 GeV, At ^ 0.0 GeV^ with an 
RMS of 0.02 GeV^. 

The final set of filters uses kinematic cuts to reduce the number of background events. 
Each event was required to have: 

• a total energy in the forward, inner part of the LAr calorimeter (1.7 < rj < 3.3), of 
less than 400 MeV suppressing events in which the proton dissociates. 
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Figure 1: (a) MC reconstructed versus generated W^p] (b) MC reconstructed versus generated 
t for the low W-yp analysis. Only elastic events are shown. The distributions are shown after 
all cuts were applied. 

• pt > 0.5 GeV for the high pt track and pt < 0.4 GeV for the low pt track, matching 
the region of flat efficiency in pt for the DCR0 trigger, 

• I tree I < 0.5 GeV^, with t^-cc = —Pt'^ being the reconstructed t, reducing the number 
of dissociative and inelastic events, 

• W^p'^ > 40 GeV, cutting most of the proton beam-gas events present in data and 
restraining the measurement to the range where there is reasonable acceptance, 

• Mtttt e [0.52, 1.17] GeV, hmiting the measurement to the range where there is a 
reasonable signal to background ratio. 

A total of 358 coUiding bunch events and no pilot events survive all cuts. 

The trigger efficiency was measured using independent samples of events, where there 
was a chance overlap between the reaction ep ep^p and the Bethe-Heitler process 
ep e'-fp. These events were triggered by the detection of the electron from the Bethe- 
Heitler process in accidental coincidence with a signal from other detector components. 
The LAr veto was found to be 100% efficient for the selected events. Triggers independent 
of the CTD were used to determine the vertex efficiency e^-vertex- With shghtly relaxed 
selection criteria, e^-vertex = 0.72 ± 0.03(stat.) was obtained, in agreement with MC 
simulations. The relative efficiency of the DCR0 and remaining requirements for events 
already fulfilling the 2;-vertex trigger, encR.^,, was determined from a sample triggered 
by the electron and the z-vertex condition only. Applying the final selection procedure 
and cuts ensured the same event properties as for the signal sample and gave eocR^ = 
0.29 ± 0.06(stat.). The results were multiplied together to obtain the overall trigger 
efficiency of e^ng — 0.21 ± 0.04. The error quoted is dominated by the low statistics in the 
overlap sample. 
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Backgrounds in the final sample were estimated using the PYTHIA MC generator. 
Only the diffractive channels were found to contribute significantly. The estimated back- 
ground fractions in the final sample of p*^ candidates was 13 ± 7 % from GD, 9 ± 5 % from 
PD, 10 ± 5 % from DD, ± 1 % from ND and 1 ± 1 % from elastic production of uo and 
meson^. The errors on the background includes the uncertainty in the decomposition of 
the total photoproduction cross section. 

Figure || shows the dipion mass spectra of the accepted data (without background 
subtraction) for the low W^p analysis along with the MC prediction. The MC mass 
shape, containing a Ross-Stodolsky skewed mass distribution for the elastic events (see 
section |6.1| for details), matches the shape of the data histogram well. 
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Figure 2: Reconstructed M,r7r distributions for the low W^p data sample before background 
subtraction and acceptance correction. The open histogram shows the MC (PYTHIA) predic- 
tion for all possible channels (normalised to the data). The shaded histogram shows the MC 
prediction for all channels apart from EL. 



5 The High W^p Analysis 

The data for this measurement were collected during a short dedicated run in the 1994 
data taking period with a vertex shifted forward of its nominal position by ~ 70 cm. With 
this vertex position the acceptance for particles traveUing backwards at small angles to the 
beam pipe is enhanced. The HERA machine was operated with 153 colliding bunches of 

^The pt cuts imposed on the tracks strongly suppress events from the reaction 7p — s- (f>p, 
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27.5 GeV positrons and 820 GeV protons with in addition 17 proton and 15 positron pilot 
bunches. The data taken correspond to a total integrated luminosity of 23.8 ± 1.3 nb~^. 

Photoproduction events were triggered by the coincidence of a signal in the electron 
tagger {E'^ > 4 GeV) and a signal in the ToF system coming from the time interval 
expected for ep interactions. This trigger has been shown to have a high (35 %) and well 
understood acceptance for all photoproduction processes for data taken over a similar 
W^p range §. 

The main source of background in the triggered sample was found to be electron beam 
interactions with residual gas or with material inside the beam pipe. This background 
was reduced by requiring at least one BPC hit and a total energy deposition greater than 
0.2 GeV in the BEMC or LAr calorimeters. The measurement was further restricted to 
the kinematic range 0.3 < y < 0.5 where there was good electron tagger acceptance. 

Two classes of diffractive events were distinguished on the basis of a 'rapidity gap' i.e. 
a region of laboratory pseudorapidity where no particles are observed in the final state. 
The first, termed an elastic proton sample where in most cases the proton remained intact 
(GD and EL), was selected as follows. No activity was required above threshold in the 
forward detectors (FMD and PRT). The pseudorapidity, r^max, of the most forward energy 
deposit greater than 600 MeV in the LAr or BEMC calorimeters or the most forward track 
in the CTD was required to be less than 3. The second class of event, where the proton 
dissociates (DD and PD), was selected by requiring there to be signals in either FMD, 
PRT or r^max > 3. ND events in the proton dissociation sample were suppressed by 
requiring Arj, the largest region in the detector where no tracks or energy deposits were 
found, to be greater than 2 units of pseudorapidity with the upper edge of the rapidity 
gap having r] > 2.5. It should be noted that these latter requirements also suppress those 
events where the proton fragments into a high mass system and in which the rapidity 
gap is small or non-existent. The distinction between these DD events and ND events is 
experimentally not well defined. Studies with the MC generator indicate that these cuts 
restrict the proton dissociation sample to a proton dissociation mass of My < 10 GeV. 

The four vector {Eh, Pxh, Pyh, Pzh) of the hadronic final state excluding the proton 
or the dissociative proton system was determined by combining calorimeter (LAr and 
BEMC) and tracking (CTD) information in a procedure that avoided double counting and 
optimised the resolution. For the proton dissociation sample, only those tracks/energy 
deposits at angles more backward of the rapidity gap were included. The reconstructed 
invariant mass of the hadronic final state excluding the proton, Mj?'^, was then calculated 
making the assumption that P|/j + P^^ ~ |t| -C {Mx^Y and using the approximation 



Here E^ is the electron beam momentum and y is reconstructed from the final state 
electron. This method resulted in good resolution, as demonstrated by figure |^, even 
though a significant proportion of the hadronic energy is not detected in the backward 
region. It should be noted, however, that the resolution width in is much greater 
than the width of the p° decay, making it impossible to determine the line shape with 
this analysis. 



y = {Eh- P,h)/{'2E,) so that 




(4) 
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Figure 3: MC reconstructed versus generated Mx for the high W^p analysis. Ah diffractive 
channels (EL, GD, DD, GD) are shown. 



A sample of events where elastic p° production dominates was selected by requiring 
0.4 < Mx < 1.26 GeV in the elastic proton sample. 

Background in the elastic p° sample was observed and attributed to four dominant 
sources: electron beam induced, GD, PD where none of the proton fragments were de- 
tected and EL from u and production. 

The electron beam induced background was estimated using pilot bunch events, was 
found to contribute ~ 20 ± 5 % and was removed from the sample. 

The MC generators (PYTHIA and PHOJET) were used to estimate the fraction of 
GD, PD, EL, DD and ND background. The relative cross sections used in the generation 
for the four diffractive processes described above were chosen by comparison of the 
spectra in data and MC for both the elastic proton and proton dissociation samples. The 
region where > 1.26 GeV was found to be dominated by GD and DD and that 

where M^'^ < 1.26 GeV by EL and PD. Since the relative cross sections obtained in this 
comparison for PYTHIA and PHOJET were in broad agreement, the same fractions were 
used for both MC generators (see equation This choice also satisfied the constrains 
imposed in |Q . Figures ^ (a,b) show comparisons of the data with MC predictions (with 
the cross section ratios as described above) for the M^^ distributions of the elastic proton 
and dissociated proton samples. For these plots and subsequent results the average of the 
two MC generators was used. 

The ratio of GD events produced in the MC generators for 0.4 < Mx'^^ < 1.26 GeV to 
EL events was found to be ~ 6 %. This number agrees well with measurements made at 



lower energies |£0|- The GD, PD, EL, DD and ND background fractions in the elastic p 







sample were estimated to be 12 ± 6 %, 6 ± 3 %, 10 ± 5 %, 3 ± 2 % and ± 1 % respectively. 
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Figure 4: Reconstructed Mx distribution for the high W^p data for a) the elastic proton selec- 
tion and b) for the proton dissociation selection. The solid histograms show the MC prediction 
for all possible channels. The dashed histograms show the MC prediction for a) all channels 
apart from EL and b) all channels apart from PD. The MC distributions were normalised to 
the data in the region Mx < 1-26 GeV in the elastic proton selection. 



As for the low W^p case the errors on the background estimates accommodate the range 
of uncertainty in the decomposition of the total photoproduction cross section into its five 
constituent processes 0. All backgrounds were subtracted statistically from the data. 

The acceptance for the electron tagger and for the detection of the hadronic final 
state were determined separately. The first was determined in a procedure identical to 
that described in 0, was found to be 61%. The hadronic final state acceptance was 
determined by taking the mean of the values obtained from two MC generators, namely 
39 % and 40 % for PYTHIA and PHO JET respectively. 



6 Results 

6.1 The Differential n'^7r~ Cross Section dcr/dM^Tr 

The corrected differential cross section of the invariant mass of the tt+tt^ system, dcr/dM^rTr, 
was determined for the data taken at (W^p) = 55 GeV from the raw distribution (shown 
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in figure 0) using 

da _ Nun - iVbgd 1 1 



dM,, $7££trig£bmC l + i?PD AM 

Here: 

- A^bin is the number of events reconstructed in eacli mass bin, 

- iVbgd is tlie estimated number of background events in the bin, from GD, DD, ND 
and EL, 

- RpB is the estimated ratio of PD to EL (taken as a constant scahng factor for all 
bins), 

- $/ is the photon flux found by integration of equation ||, for 40 < W^p < 80 GeV 
and gLn<Q'<0.5 GeV^, 

- C is the integrated luminosity, 

- £trig is the global trigger efficiency as calculated in section ^, 

- £bm is the bin by bin selection efficiency given by the MC before the trigger selection, 
averaged over the kinematical range generated, 

- C is the efficiency correction equal to the ratio of data to MC track reconstruction 
efficiencies, 

- AM is the bin width. 

The acceptance, eun, varied between 2 % and 37 % with an average value of 15 %. 
The integrated photon flux is = 0.0621. The production cross section is shown in 
flgure ^ The mass dependence of dcr/dM^7r is consistent with a large contribution from 

photoproduction. The peak is shifted to lower values than the nominal p° mass and 
there is an enhancement (suppression) of the cross section at values lower (higher) than 
the nominal p° mass. This asymmetric shape, or skew, of the distribution is a well known 
feature of p° photoproduction which can be attributed to an interference between the 
resonant and non-resonant production of two pions. In the phenomenological approach 
by Ross and Stodolsky the differential cross section is assumed to follow 

'^'^ -- f,-BW,{M^^)-{M,/M^^Y^^ (6) 



dM 



TTTT 



where /„ and are constants and 



where Mp is the p° mass and where the momentum dependent width, following the sug- 

taken to be 



gestion of Jackson p2|, is taken to be 



Here Fq is the p° width, q* is the vr momentum in the tt+tt^ rest frame and is the 
value of q* when M-j^t^ = Mp. The distortion of the mass spectrum is characterised by the 
Ross-Stodolsky parameter riRs, which becomes zero for production of just the resonance. 
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Figure 5: The da/dMn-jr distribution at (Wjp) = 55 GeV for \t\ < 0.5 GeV^. Only statistical 
errors are shown. The solid curve represents a fit to the distribution assuming the parameterisa- 
tion of Ross and Stodolsky (equations 7 and 9) . The dashed curve shows the pure Breit Wigner 
(equation 7 with n^s = 0). 



The spectrum was fitted, using the parameterisation of equation over the range 
0.52 < Mtttt < 1.17 GeV with quantities fp, Mp, Fq and n^s left as free parameters. The 
fit is shown in figure ^ and the results are summarised in table |l|. The values of Mp and 
Fq are in agreement with particle data group (PDG) values of Mp = 769.9 ± 0.8 MeV and 
Fq = 151.2 ± 1.2 MeV [^. The value of the differential cross section at the nominal p° 
mass was measured to be 



da 



dM^ 



44 ± 5 (stat.) ± 12 (syst.) /ib GeV" 



--M„ 



The systematic error was determined in a similar way to that described in section 
Spital and Yennie have suggested that the value of the differential cross section at the 
nominal p° mass offers a way to determine the total elastic p° cross section independently 
of the assumptions made for the form of the line shape since it is at the nominal p° mass 
where there is little or no interference of any additional coherent production [^. With 
their assumption of a standard Breit Wigner for the p° resonance line shape, we get 



vrFo dcr 



10.6 ± l.l(stat.) ± 3.0(syst.) fih 
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Parameter 


Value 


M, 


783 


± 13 


MeV 


To 


153 


± 19 


MeV 


/p 


6.23 


± 0.62 




nns 


5.84 


± 1.05 






12. 


81/12 





Table 1: Results of the fit of the da/dM^^ distribution for (W^p) = 55 GeV and \t\ < 0.5 GeV^ 
assuming a Ross-Stodolsky skewing factor (see text). The errors shown are statisticah 



In order to investigate the sensitivity of the results to changes in the assumed pa- 
rameterisation of the mass spectrum, the data were reanalysed with a fit based on the 
Soding interference model |2^ and with different assumptions for the momentum depen- 
dent width. In the Soding model the skewing of the mass spectrum is explained by the 
interference of a resonant —>■ tt+tt^ amplitude and a p-wave tttt Drell type background 

p7[ can be written as 



term ||26|. A simplified parameterisation of the model used in 



da 



fp-BWp{M^^) + fi-I{M^ 



(9) 



where BW plM^^j^) is defined in equation |^ and the interference term is 



Ml 



(10) 



The quantities /p, /j, Mp and Vq are left as free parameters. The fit was also done for 
both models with the following two parameterisations of the momentum dependent width, 
used by some other experiments |2^ 



Tp{M^^) = 



) = ro(^f, 
% 



(11) 

(12) 



The results of the fits are summarised in table 0. All of the fits yielded values for Mp and 



Fq in agreement apart from that using the Soding model with the width of equation 11 



6.2 The Elastic Photoproduction Cross Section 

The 7p n^n^p cross section at {W-yp) = 55 GeV was obtained by integrating the fitted 
form of the differential invariant mass spectrum of equation ^ with a momentum dependent 
width as in equation |^. This procedure involves an extrapolation outside the measured 
region, 0.52 < M^^ < 1.17 GeV, to the region 2M^ < M^^^ < Mcut, where, following the 
procedure used by 0, the upper mass limit is taken to be Mcut = Mp + 5To. In making 
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Model 


/ndf 


Mp /MeV 


To /MeV 


f„ /ub 


' ^Ho 


fr //ibGeV-i 


RS +(8) 


1.07 


783 ± 13 


153 + 19 


6.23 + 0.62 


5.84 + 1.05 




RS +(11) 


1.09 


791 ± 14 


163 + 23 


6.11 + 0.65 


6.84+1.00 




RS +(12) 


1.09 


784 ± 13 


153 + 19 


6.20 + 0.62 


5.95 + 1.01 




S +(8) 


0.84 


787 ± 12 


176 + 28 


6.70 + 0.55 




3.39 + 0.85 


S +(11) 


1.02 


815 + 16 


205 + 40 


5.86 + 0.49 




5.07+1.20 


S +(12) 


0.88 


793 ± 13 


181 + 30 


6.57 + 0.54 




3.74 + 0.92 



Table 2: The results of the fits of the da/dM-j^.,^ spectrum for the low W^yp analysis. In the 
column marked 'Model': RS is the Ross-Stodolsky parameterisation of equation 6, S is the 
Soding parameterisation of equation 9. The bracketed numbers refer to the equation describing 
the momentum dependent width. The errors quoted are only statistical. The column labels are 
specified in the text. 



the extrapolation we assume a t dependence with a slope parameter b = 10.9 GeV ^, as 
obtained from the fit of the dcr/dt distribution in section |6.3| , so that 



dM^^ dt 



A{M^^)-e''\ (13) 



The integration is performed over the limits —0.5 GeV^ < t < tmin, where |tmin| is the 
minimum value of \t\ which is kinematically allowed. The fraction of the cross section 
that lies outside the measured kinematic region in M^^ , ^, is found to be ^ = 0.15. The 
resulting total 7p — ^ tt^tt^p cross section is 

cr(7p TT+TT-p) = 11.2 + 1.1 (stat.) + 3.1 (syst.) /ib. 

The resonant part of the total cross section is obtained by integrating the function 
fp ■ BW p{Mt,t,) (i.e. taking urs = in equation ^), assuming the same t dependence 
as before, over the prescribed range in M^^. The value of the cross section, measured at 
{W^p) = 55 GeV for 2M^ < M^^ < Mp + 5To and -0.5 GeV^ < t < t,nin is then 

a{^p^pOp) = 9.1 + 0.9(stat.) +2.5(syst.) /ib. 

At the present level of experimental precision there are no corrections to be made for 
p° decay modes other than to vr+vr". 

For the chosen value of M^nt = Mp + 5To the extracted cross section, (t(7p P^p), was 
found to be very stable against changes in the parameterisation of the da/dM^^T^ spectrum. 
Allowing riRs to vary within the statistical errors resulted in a change of ±6.0 %. Repeating 
the fit with values for Mp and Fq fixed to the PDG values gave an increase of 7.0 %. Using 
the parameterisation of equation |^ gave an increase of the cross section of 5.9 % and the 
maximum variation caused by different assumptions of the momentum dependence of the 
width was 6.6 %. 

If, however, we do not impose an upper Mtttt limit but use the parameterisations to 
extrapolate over the full kinematic range available (i.e. Mcnt = W.yp — Mp) as was the 
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procedure for measurements made at lower energies, we find a very large variation in the 
cross section extracted with the various methods (see table ^ for a summary of all the fits). 
This is because some of the parameterisations assume the p*^ resonance to have a long 
tail extending up to high values of M^^. Unlike the situation at lower energies kinematic 
constraints do not suppress this tail to levels where it can be neglected. 



Model 




M,ut = W^p - Mp 




O-TTTT / /ib 


apljjh 






RS +(8) 


11.2 


9.1 


11.2 


10.5 


RS +(11) 


11.4 


9.3 


11.4 


15.9 


RS +(12) 


11.5 


9.2 


11.5 


13.5 


S+(8) 


11.2 


9.6 


10.8 


11.4 


S +(11) 


12.3 


8.5 


17.0 


14.5 


S +(12) 


11.5 


9.6 


14.5 


14.8 



Table 3: The extracted o"(7p 'k~^it~p) and (t(7p p^p) cross sections from the fits of the 
dc7/dM,r7r spectrum for two choices of Mcut (see text) for the low W-yp analysis. In the col- 
umn marked 'Model' RS is the Ross-Stodolsky parameterisation of equation 6, S is the Soding 
parameterisation of equation 9. The bracketed numbers refer to the equation describing the 
momentum dependent width. 

Since it was not possible to determine the Mtttt and t dependence of the cross section 
at {W^p) = 187 GeV, we took the values of n^s and b found in the low W^p analysis^] 
to reweight the MC distributions used for the acceptance calculation and to determine 
numerically ^ (the fraction of the cross section outside the measured region) and the 
ratio of the cross sections rnr = cr(7P P^p)/o'{lP n'^ii'p). The values obtained were 
= 0.05 and r^r = 0.82. The elastic cross section was then determined using 

a{^p ^ p^p) = ^^^^ . __ . r^^^ (14) 

where N is the measured number of events in the mass range 0.4 < Mtttt < 1.26 GeV, A^bgd 
is the estimated number of background events, $/ = 0.00903 the photon flux found using 
equation ^ £ is the integrated luminosity and e is the acceptance. The elastic jp —>■ p^p 
cross section for -0.5 GeV^ < t < tmin and 2M^ < M^^ < Mp + 5To at (W^p) = 187 GeV 
is then 

(y{lP^P^p) = 13.6 ± 0.8 (stat.) ±2.4 (syst.) /ib. 

The breakdown of the various contributions to the systematic error of the elastic p'^ 
cross section measurement at (W^p) = 55 GeV for the mass range 2M^ < M^rTr < Mp + 5ro 
can be found in table |[ The primary source of systematic error comes from the uncertainty 
of the trigger efficiency, where the large error of 20 % is due to the low statistics in 

■^This procedure is likely to be valid since the measured value of nRs is close to values obtained at 
lower energies j27[ |8[ |o), indicating there is little dependence of nns on W-yp. 
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the overlap sample of events used to determine it. An error of 5 % is assigned to the 
uncertainty in the scanning yield, which was estimated by taking the spread in the yield 
as found by several scanners. The track fitting error of 4 % was estimated by taking the 
difference between the acceptances found in data and MC. An error of 3 %, arising mainly 
from limited statistics, was assigned to the MC acceptance. Changing the t dependence 
within the errors of the measurement gave a 12 % error. To estimate the error due to 
background contributions each background was varied by ±50 % and the M^^^^ distribution 
re-fitted. An error of 6 % was assigned to reflect the spread in the extracted cross section 
due to different assumptions for the parameterisation of the M^^r dependence. The error 
in the integrated luminosity measurement was 5 %. 



Source 


Error 


Trigger efficiency 


20% 


Scanning yield 


5% 


Track fit efficiency 


4% 


MC acceptance 


3% 


CD background 


8% 


PD background 


5% 


DD background 


2% 


ND background 


1% 


EL background 


1% 


t dependence 


12% 


Resonance extraction 


6% 


Luminosity 


5% 


Total systematic error 


28% 



Table 4: Sources of systematic errors for the elastic cross section measurement at (Wyp) = 
55 GeV. 

For the high W.yp analysis the biggest uncertainty on the elastic cross section mea- 
surement comes from background in the elastic p° sample. A systematic error was taken 
of 50 % of each background estimated using the MC simulations and 25 % for the electron 
beam gas background. The error on the electron tagger efficiency was 5 % . An error 
of 6 % was assigned for the MC acceptance being due to the difference between PHOJET 
and PYTHIA and the MC statistics. An error of 1 % was found by varying the t slope 
parameter by ±4 GeV~^ around the central value of 11 GeV^^. The 20 % overall hadronic 
energy scale uncertainty of the BEMC pl| gave rise to a 4 % error. A 6 % error was as- 



signed to account for the uncertainty in the method of extraction of the resonant cross 
section. This error was estimated by taking the spread in the value of the cross section 
obtained using the results of each fit to the low W.yp data as described above, reweighting 
the input MC M^^r distributions and recalculating the acceptances and values for ^ and 
Tnr. The error in the integrated luminosity measurement was 5 %. A breakdown of the 
systematic error is shown in table 
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Source 


Error 


(tT) nfi CTTon ri n 


O /(J 


PT^ n^iPKOTonrin 


4% 


T^T^ n^iPKOTonnn 


2% 


IMT) nfiPK o"rrniTin 


1 % 




7% 


MC acceptance 


6% 


BEMC hadronic scale 


4% 


t dependence 


1% 


e-gas background 


5% 


e-tagger acceptance 


5% 


Resonance extraction 


6% 


Luminosity 


5% 


Total systematic error 


18% 



Table 5: Sources of systematic errors for the elastic (P cross section measurement at (W^p) = 
187 GeV. 




Figure 6: The dependence of the elastic photoproduction cross section, 17(7^ p^V)-, on 
VF-yp for HI and previous measurements 1^, |6|, |3^, The inner vertical bars on the HI points 
denote only the statistical error, the outer ones contain statistical and systematic errors added in 
quadrature. The horizontal error bars show the range in W-^p over which the measurements were 
made. The solid curve shown is a curve based on a soft pomeron calculation, the normalisation 
of which is fixed by the data at low energy [ 11 1 . 
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Figure |^ summarises the dependence on Wjp of measurements of the elastic photo- 
production cross section 0, |, ^ |3^ including our new measurements. Also shown is the 
theoretical prediction based on the VMD model and non-perturbative Regge theory [|ll[ 
(solid curve). The new high energy measurements at HERA demonstrate that the W^p 
dependence of the cross section is consistent with the slight increase characteristic of this 
model in which the leading exchange is a 'soft' pomeron. 



6.3 The t Dependence of the Elastic p Photoproduction Cross 
Section 

Diffractive scattering is characterised by a low momentum transfer between the vector 
meson and the proton, resulting in a steeply falling distribution in t. This is quantified by 
the slope parameter, b, determined by fitting the form e''* to the t dependence of da/dt 
for the data taken at (W^p) = 55 GeV. 




- b = 10.9 ±2.4 [GeV"^] ^ ^ . 

d0/dt(O) = 79.9 ± 30.5 [^ib/GeV^] 

'"' I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

-t [GeV^] 

Figure 7: The differential cross section, da/dt, at {W^p) = 55 GeV is fitted over the range 
0.025 < \t\ < 0.25 GeV^ with the function da/dt = ae^^ (only statistical errors are shown). 

The cross section da/dt was determined by integrating the mass spectrum over the 
measured range (0.52 < M^tt < 1-17 GeV) and correcting for losses and backgrounds using 
the MC simulation. A correction was made for losses in the tails making the assumption 
that nRs, the Ross-Stodolsky parameter, remained constant with t, so that ^ = 0.15 at 
each value of t. The results are shown in figure |^. The data for 0.025 < \t\ < 0.25 GeV^ 
are well parameterised with a simple exponential with slope parameter 

b = 10.9 ± 2.4 (stat.) ± 1.1 (syst.) GeV"^ 
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and intercept 



dt 



79.9 ± 30.5 (stat.) ± 21.9 (syst.) ^b/GeV^ 



The systematic error arises from the variation of the MC input parameters for the skewing 
of the hne shape and the h slope within the statistical errors of the fits and also by 
allowing each background to vary by 50 %. For the intercept the systematic errors on 
the normalisation are also taken into account (see section |6]^). For \t\ > 0.25 GeV^, 
deviations from a simple exponential dependence were shown by other measurements of 
diffractive vector meson photoproduction ^ 

The measurement of h is compared to previous determinations in figure |. The recent 
ZEUS measurement |^ is shown along with lower energy data as compiled by Aston et 
al. 113, |8|, |9|, |0|, |33. The ZEUS and HI results together show that the shrinkage of 
the t dependence of da /dt continues into the HERA energy range as expected in the soft 
pomeron picture. 
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Figure 8: Determinations of the exponential slope parameter, 6, as a function of E^^ the 
photon energy in the rest frame of the proton. The inner vertical bar on the HI point denotes the 
statistical error only, the outer one contains statistical and systematic errors added in quadrature. 
The horizontal error bars show the range in over which the measurements were made. The 
curve is a prediction based on pomeron exchange ||TT|. 



6.4 Polar Angular Decay Distribution of the 

A feature of the diffractive production of hadrons, in addition to a slight dependence of 
the cross section on W^p and a peripheral t dependence, is s-channel helicity conservation. 
This is best studied in the helicity frame, defined as the rho meson rest frame, with the 



quantisation axis along the direction of the rho in the 7p frame |2^. In this frame, the 
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Figure 9: The acceptance corrected cos 9* density distribution at (W^p) = 55 GeV. The data 
are fitted to the function in equation 16. The errors quoted are only statistical. 



polar angle distribution of the tt"'" according to the formalism in |^ is 

oc(3r°o^-l)cos2r + (l-0, 



d cos 9* 



(15) 



where Tqq is the spin density matrix element which specifies the probability that the p° 
meson is produced in the spin substate 0. 

In order to determine Tqq, the distribution in cos0* was obtained using the same selec- 



tion, subtraction and correction procedures as for the da/dt analysis. The expression |15 
was then fitted to this distribution to obtain a value of 



,04 
00 



-0.11 ± 0.11 (stat.) ± 0.04 (syst.) 



Figure ^ shows dcr/dcos0* and the fit of the expression |T5|. In the helicity frame the 
p° is thus observed to be produced predominantly in spin substates ±1 and not in spin 
substate 0, consistent with the hypothesis of s-channel helicity conservation at = 0. 



7 Conclusions 

With the HI detector we have measured the cross section for the elastic photoproduction 
of p° mesons and found cr(7p p^p) = 9.1 ± 0.9 (stat.) ± 2.5 (syst.) /ib at (W^p) = 
55 GeV and a{-fp p^p) = 13.6 ± 0.8 (stat.) ± 2.4 (syst.) ph at {W^p) = 187 GeV for 
2Mt, < Mtttt < Mp + 5Tq. At the former W^p, the slope parameter of the distribution in t is 
10.9±2.4 (stat.) ± 1.1 (syst.) GeV~^ and the decay polar angular distribution is consistent 
with a pure sin^ 9* distribution. We thus verify the extension to HERA energies of the 
properties found at lower energies for elastic photoproduction of the meson, namely a 
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steep forward distribution of the produced p° mesons with a continued shrinkage of the t 
dependence with increasing W^p and a decay angular distribution that is consistent with 
s-channel hehcity conservation. A slow rise of the elastic p° cross section with Wjp, as 
expected in a picture of interactions based upon the soft pomeron which also matches the 
rise of the total photoproduction cross section, is shown to be consistent with observation. 
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